Kinetic Derivation of the Hydrodynamic Equations for Capillary Fluids 



o 
o 

(N 

>>: 



S. De MartinoQ M. FalangaQ and S. I. Tzeno\@ 

Dipartimento di Fisica -Universita degli Studi di Salerno, Via S. Allende, 
Baronissi (SA), 1-84-081 Italy; INFM umtd di Salerno, 84081 Baronissi (SA), 
Italy; INFN, Sezione di Napoli, Gruppo Collegato di Salerno. 
(Dated: May 2004) 

Based on the generalized kinetic equation for the one-particle distribution function with a small 
source, the transition from the kinetic to the hydrodynamic description of many-particle systems is 
performed. The basic feature of this new technique to obtain the hydrodynamic limit is that the lat- 
ter has been partially incorporated into the kinetic equation itself. The hydrodynamic equations for 
capillary fluids are derived from the characteristic function for the local moments of the distribution 
function. The Fick's law appears as a consequence of the transformation law for the hydrodynamic 
quantities under time inversion. 



PACS numbers: 05.20.Dd, 47.10.+g, 05.40.Jc 

Keywords: capillary fluids; kinetic equation; hydrodynamic picture 



i 

q, 

O ■ 



> 

O 
in 
o 

o 

o 

Or 

> : 
X- 



I. INTRODUCTION 

The Van der Waals gradient theory [| was originally 
developed at the end of the 19-th century as an "effec- 
tive" physical picture to describe the critical region. This 
theory considers a local Helmholtz free-energy composed 
of two parts. The first one describes the homogeneous 
behaviour of the system, while the second part character- 
izes the non-homogeneous one. In addition, the first part 
is proportional to the fluid density g while the second 
one is proportional to the square of the density gradi- 
ent. Much effort has been devoted to apply the Van der 
Waals gradient theory to various fluids. At present, it 
can be considered as the most simple and comprehensive 
model in the physics of interfaces and capillarity. The 
associated hydrodynamic equations can be written as 

" ,J - ' " (1) 
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dt 
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where !F(g, a) is a function of the density g and of 
a = (l/2)|V ( o| 2 and V is the current velocity. This 
formulation of the Van der Waals theory was originally 
due to Korteweg 3] , who proposed a continuum mechan- 
ical model in which the Cauchy stress tensor apart from 
the standard Cauchy-Poisson term contains an additional 
term defined as 



(3) 



T = [-p + aV 2 g(x;t) + p\Vg(x;t)\ 2 )l 
+5Vg(x; t) <g> V£>(x; t) + 7 (V <g> V)q(x; t), 



where 1 is the unit tensor. As already mentioned by 
Dunn and Serrin 4] , the modern terminology concerning 
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the Korteweg model refers to elastic materials of grade n, 
where the particular case of n = 3 has been well studied 
in recent years 

Equations (QJ and (JSJ) have been linked recently QJ 
to a nonlinear Schrodinger equation viewed as a partic- 
ular set of hydrodynamic equations describing the so- 
called nonlinear Madelung fluid 0, 0]. This link be- 
tween capillarity and the Schrodinger equation can shed 
more light onto the so-called quantum-like approach to 
many-particle systems such as beams in particle accelera- 
tors and beam-plasma systems. The standard procedure 
in this direction is to approximate the physical systems 
characterized by an overall interaction with a suitable 
Van der Waals mean field theory. 

To avoid misunderstanding, it is worthwhile to note 
that Schrodinger equation alone does not provide an en- 
tire quantum mechanical picture. It should be necessar- 
ily supplemented by a theory of quantum measurement 
and consequently by a proper physical interpretation of 
wave packets. In the quantum-like approach, the many- 
particle systems are described in an effective way as a 
whole. Based on the above considerations, it appears 
interesting and attractive to explore the possibility of a 
rigorous derivation from kinetic theory of the general hy- 
drodynamic picture thus discussed. 



II. THE GENERAL FRAMEWORK 

The starting point of our analysis is the equation for 
the microscopic phase space density Nj^(x,p;t) 



M 



dt 



m 



( P N M ) + d p ■ [F M {x,p;t)N M ] = 0, (4) 



for a system consisting of N particles, which occupies 
volume V in the configuration space. Here x and p are 
the coordinates and the canonically conjugate momenta, 
m is the particle mass and Fm(x, p; t) is the microscopic 
force, which apart from the external force includes a part 
specifying the type of interaction between particles. Sup- 
pose that at some initial time £q the microscopic phase 
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space density is known to be ]Vmo(x, p; to). Then, the 
formal solution of equation |@J for arbitrary time t can 
be written as 



iV M (x, p;i) =<S(t;t )iVMo(x,p;to) ; 



(5) 



where S(t; to) is the evolution operator, specifying the 
Hamiltonian flow. 

The choice of the initial JVmo(x, p;io) is based on the 
knowledge of the microscopic characteristics of the sys- 
tem. Due to the extremely complex particles' dynam- 
ics, full consistent description is not feasible. Therefore, 
the detailed information on the microscopic level is in- 
complete. If our system is a complex one in the sense 
that both the external forces and the collective forces are 
highly nonlinear, a dynamic instability of motion is likely 
to occur on a characteristic time scale r. The only in- 
formation available to an outside observer by means of 
a macroscopic measuring device is a coarse-grained den- 
sity distribution with a smoothing function, which takes 
into account the dynamic instability of motion. Thus, we 
assume 

iY MO (x;to) = N M (x;t ) = [ d 3 zG(x;t | z)N M {z;t ), 



(6) 

where for simplicity the explicit dependence on the mo- 
mentum variables p has been suppressed. To take into 
account the initial preparation of the system, one has to 
displace the initial time to at — oo and perform an aver- 
age over the pas t history of the system. Then equation 
(gj becomes [IcJ 



dN 



M 



at 



-V-(pN M )+ 9 p -[F M (x,p;t)7V M ] 
m 

= ^N M -N M y (7) 



Since the collision time is supposed to be much smaller 
than the time r, the collision integral can be dropped 
and the kinetic equation for the one-particle distribution 
function /(x, p; t) can be written as 



dt 



-•V/ + F(x,p;t)."£ p / =-(/-/ 
m t K 



(8) 



The right-hand-side of equation (JSJ is regarded as a " col- 
lision integral" , and it can be represented as 



^ oo i 



(-D' 
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(9) 



where 



(10) 



with Az = z x. As a first very interesting step, we 
consider the diffusion approximation 

1(7 ~f) =-V k [A k (x;t)f] + ±V k Vi[B kl (x;t)f}, (11) 
where 

A fe (x;i) = -(Az k )<?l B kl (x;t) = -(Az k A Zl )^, 

T T 

(12) 

and a summation over repeated indices is implied. For 
Hamiltonian systems the well-known relation 



A k (x;t) = -ViB u (x;t) 



(13) 



holds, which gives A = for B k i — const. 

In passing, it is worthwhile to mention that using the 
principle of maximum information entropy formulated by 
Jaynes, it can be shown [Tfl | that the smoothing function 
G(z; t\ x) is of the form 



G(z;t|x) 



1 



r 3/2 



x exp 



-I z - ;z)„ , | C (x;t)(z - Wx _ t 



detC 

(G) 



The quantities 
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,(G) 
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{z k zi 



,(G) 
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(14) 



(15) 



are the first and the second moment of z at the instant 
of time t + r, provided that z measured at the instant 
t equals x [i.e. z(t) = x]. In addition, C(x;t) is the 
covariance matrix defined as 



Cfcz(x,t) = 2 



MS 



(^)i G * ) (^>S ) 



(16) 



The generalized kinetic equation JSJ has a form anal- 
ogous to the Bhatnagar-Gross-Krook (BGK) equation, 
widely used in the kinetic theory of gases hll. There 
is however, an important conceptual difference between 
the two equations. In the BGK equation the function / 
should be replaced by the equilibrium distribution func- 
tion fo describing the global equilibrium and the charac- 
teristic time t should be replaced by the corresponding 
relaxation time. The smoothed distribution function in 
equation (JHJ characterizes a local quasi-equilibrium state 
within the smallest unit cell of continuous medium, while 
t is the corresponding time scale. 



III. THE HYDRODYNAMIC APPROXIMATION 

Rather than following the standard approach in deriv- 
ing the hydrodynamic picture, we introduce the charac- 
teristic function 



S(x,w;t) = J d 3 p/(x,p;i)e-™"P, 



(17) 
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instead. It is straightforward to verify that Q satisfies the 
following equation 

^ + -\7-l) w g + iwFg = -V-(A£) + iv„V s (B„ s £). 
ot m 2 

(18) 

Note that the local moments (p™ 1 ^ 2 ■ ■ ■P k k ) can be ob- 
tained from the characteristic function according to the 
relation 



(19) 



where n± + Ti2 + • • • + n k = I- The well-known hydrody- 
namic quantities, such as the mass density g, the mean 
velocity V( + ) of a fluid element and the hydrodynamic 
stress tensor life; can be defined as 



£>(x; t) = mnQ (x, 0; t) = ran 



yd 3 p/(x,p;i), (20) 



e(x;t)V (+) (x;t) = ind^Q 



w=0 



n / d 3 pp/(x,p;t), 
(21) 



n fci (x;i) 



n a 2 g 



ra dwkdwi 



n 
m 



J d 3 pp kP if{x,p;t), 
(22) 

Here, n = lim (N/V), implies the thermodynamic 

N,V— >oc 

limit . Defining also the deviation from the mean ve- 
locity as 



mc (+) = P mV (+h 
and using the evident relation 



(23) 



J d 3 pc (+) (x,p;f)/(x,p;f) = 0, 

we can represent the stress tensor H mn according to the 
relation 



II mn = £>V"(_|_) m V(_|_)„ + Vmn- 



(24) 



Here 



V k l(x;t) =mn J d 3 pc (+)fe c (+); /(x,p;t), (25) 

is the internal stress tensor. 

Equation 1)18(1 and the one obtained after differenti- 
ating with respect to w k evaluated at w = 0, yield the 
Smoluchowski equation and the equation for the momen- 
tum balance, respectively. These can be written in the 
form 

|f + V ■ [f?(V (+) + A)] = iVfeV^e), (26) 



^{ Q v {+)k )+Vi{ e v {+)h v {+)l ) = ^F k 

-V • (AgV {+)k ) - VtVki + \vN n {B lne V {+)k ). (27) 

Let us consider the time inversion transformation speci- 
fied by |10Lll2| t — ► t = —t, x — > x = x and p — > p = —p. 
We argue that there exists a backward velocity V7_) (x, t) 
such that 



V (+) (x ! -t) = -V ( _ ) (x,t) 



(28) 



The transformed Smoluchowski equation i|26|l can be rep- 
resented according to 

|f - V ■ 0(-V { _) + A)] = -ivjkVKBwe). (29) 

Summing up and subtracting equations and (|29")l , we 
obtain the continuity equation 



§ + v. ( ,v) = o, 



and the Fick's law 



U k = -A k + —Vi{B k ig). 
2q 



(30) 



(31) 



Here 



V = 5(V (+) +V ( _ ) ), U=i(V (+) -V ( _)), (32) 

are the current and the osmotic velocity, respectively. It 
is worthwhile to mention that since the mean velocity of 
a fluid element is a generic function of time t, it can be 
split into odd and even part. Note that from equation 
it follows that V (+) = V + U, where V is the odd 
part, while U is the even part. 

Equation (|27|l for the balance of momentum can be 
written alternatively as 

dV (+)k , v XJV _F k 

7T, V V{-)l^lV( + ) k - 

ot y ' K ' ra 

+A l V l V {+)k - -W{P U + ^ViV n V i+)k . (33) 

Q i 

After performing a time inversion in equation i|33[) . we 
obtain 

Ot y ' v ' m 

-AtViV^k - -ViV kl - %V ; V„V(_ )fc , (34) 

Q i 

where V k i denotes the transformed internal stress tensor 
after performing the time inversion. Summing up the last 
two equations, we arrive at the sought-for equation for 
the current velocity 

dVk Fk 

+ = — + A ^iU k 

Ot m 

--yiVi k + UfViU k + ^ViV n U k , (35) 
q 2 
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where 

Vkn = \(Vkn+Vkn). (36) 

In order to find the explicit form of the internal stress 
tensor <|25[) . we observe that the maximum entropy of 
the system is realized, provided the small source in the 
generalized kinetic equation (JSJ) vanishes. This condition 
is equivalent to the condition of detailed balance in the 
case, where the collision integral (small source) is ap- 
proximated by a Fokker-Planck operator. The condition 
of detailed balance implies that the distribution function 
factorizes in the form 

/ eg ( X ,p;i) = ^)jF(p;i), (37) 

mn 

where ^"(p; t) is a normalizable function. From the above 
considerations, it follows directly that 

7> fci (x;i) = ^^(x;i)5 w , (38) 
m 

where ks is the Boltzmann constant and T is the tem- 
perature. 

In the simplest case, where the external force vanishes 
and the diffusion tensor is diagonal and isotropic, B)~i = 
I35u , we obtain 

where a = 3kgT '/m. Thus, the hydrodynamic equations 
describing a free capillary fluid have been recovered. 
In the case, where an external force is applied, the Ko- 
rteweg stress tensor contains an additional term propor- 
tional to the drift coefficient A. On the other hand from 
the principle of detailed balance, it follows that the drift 
coefficient is proportional to the external force. The 
physical implication of the latter is that the additional 
term in the Korteweg stress tensor can be regarded as a 
coupling between the external field and the mean field of 
purely hydrodynamical origin. 



IV. CONCLUSION 



Since a detailed information about the system on the 
microscopic level is incomplete, one possible way to take 
into account its initial preparation, i. e. an eventual 
dynamic instability of motion that might have set in 
and/or other large-scale characteristics, is to introduce 
a suitable smoothing procedure. As a result, the kinetic 
equation providing an unified kinetic, hydrodynamic and 
diffusion description contains a small source and is there- 
fore irreversible. Although the effective collision integral 
(small source) can be represented as a Kramers-Moyal 
expansion, for the purposes of the present paper it suf- 
fices to consider the right-hand-side of the generalized ki- 
netic equation as approximated with a properly defined 
Fokker-Planck operator. The latter form of the collision 
term is adopted as a starting point in the derivation of 
the hydrodynamic equations for capillary fluids. 
The hydrodynamic approximation is further obtained in 
a standard manner from the characteristic function for 
the local moments of the distribution function. An im- 
portant feature of the approach is that the Fick's law 
emerges naturally from the transformation properties of 
the hydrodynamic quantities under time inversion. The 
osmotic velocity is uniquely specified by the first two in- 
finitesimal moments of the smoothing function and in a 
sense is a measure of the irreversibility. 
The main result of the analysis performed in this paper, 
the hydrodynamic equations for free capillary fluids have 
been derived from kinetic theory. If an external force 
is present, the Korteweg stress tensor has to be modified 
accordingly. An additional term proportional to the drift 
coefficient emerges implying a coupling between the ex- 
ternal field and the mean field of purely hydrodynamical 
origin. 
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